Electrostatics of straight and bent nanotubes 
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Response of a single- walled carbon nanotube to external electric field, F, is calculated analytically 
within the classical electrostatics. Field-induced charge density distribution is approximately linear 
along the axis of a metallic nanotube and depends rather weakly, as \n{h/r), on the nanotube length, 
h, (here r is the nanotube radius). In a semiconducting nanotube with a gap. Eg, charge separation 
occurs as F exceeds the threshold value Fth = Eg /eh. For F > Fth, positively and negatively 
charged regions at the ends of nanotube are separated by a neutral strip in the middle. Properties 
of this neutral strip, length and induced charge distribution near the ends, are studied in detail. We 
also consider a bent nanotube and demonstrate that the number of neutral strips can be one or two 
depending on the direction of F. 

PACS numbers: 71.20.Tx, 73.63.Fg, 78.67.Ch 



INTRODUCTION 

There are two types of devices based on carbon nan- 
otubes (NTs) that are currently intensely investigated. 
These arc field-cfFcct transistors and field emitters for flat 
panel displays and x-ray sources. Although both types of 
devices were reported several years ago [l|, |2| and signif- 
icant improvement in their characterisitcs had been re- 
cently achieved (see e.g., Refs. |^^|^ and Refs. 
the understanding of performance of these devices is far 
from complete. The prime reason for this is ID-like ge- 
ometry of NT-based transistors and field emitters. Due 
to this geometry, their electrostatics is qualitatively dif- 
ferent from that in respective well-understood 2D coun- 
terparts. 

Traditionally, electrostatics for particular geometry of 
NT transistor or field emitter [l3,[l3 is studied the- 
oretically by employing a certain version of first-principle 
calculations. Notable exceptions are Refs. oaoj, m 
which classical electrostatics was used to describe the 
potential profile near the contact between metallic and 
semiconducting NTs as well as the contact between NT 
and a metal. Applicability of this description of the con- 
tact phenomena in NTs was later questioned in Ref. ■ 

In general, it is not obvious whether or not, in terms of 
electrostatics, a NT of a small radius, r, can be modeled 
by an infinitely thin sheet of electron gas wrapped into 
a cylinder. In other words, whether or not the potential 
distribution can be adequately described by the Poisson 
equation with boundary conditions imposed at distance, 
r, from the NT axis. Positive answer to this question was 
recently given in Ref. |lfil |. where the density- functional 
calculations of the extra charge distribution along the NT 
were shown to be in quantitative agreement with classical 
electrostatics analysis. 

In this situation, it is instructive to consider a model 
problem of the classical electrostatics of a NT, which al- 
lows for an analytical solution. Qualitative features of 
this solution might then yield a valuable insight into elec- 
trostatics of realistic devices. Such a problem is studied 



in the present paper. Namely, we consider a NT in exter- 
nal electric field, F, parallel to the NT axis. Separation 
of variables in the Poisson equation in this geometry is 
impossible. Still, as we demonstrate below, presence of 
a small parameter, r/h, where h is the NT length, al- 
lows one to obtain the asymptotically exact distribution 
of potential. We show that for metallic NT the density 
of induced charge changes linearly with distance from 
the NT center. For semiconducting NT with a gap. Eg, 
charge separation, which occurs as F exceeds the thresh- 
old value Eg/eh, results in formation of a neutral strip 
with a width « Eg/eF in the center of the NT. We find 
the profile of the charge density growth from the edges 
of the strip towards the NT ends. Finally, we use the 
developed approach to describe quantitatively the elec- 
trostatics of bent or wiggly NTs in external field, perti- 
nent to recent electroabsorption measurements jig , and 
demonstrate that wiggling results in multiple alternating 
positively and negatively charged regions separated by 
neutral strips. 



BASIC EQUATION 

Denote with p{z) the linear density of charge, in- 
duced by external field on the NT surface. Then the 
local value of the Fermi momentum is given by pf{z) = 
'Kh\p{z)\/2eN , where N accounts for the spin and band 
degeneracy [N = 2 or 4, and is determined by the NT 
chirality). The local chemical potential, p{z), is related 
to pf{z) via the NT energy spectrum 

fi{z) = sign(z) ^E^g/4 + v^pl{z), (1) 

where vq « 8 x 10^ cm/s is the electron velocity in 
graphene. Second relation between /i(z) and p{z) ex- 
presses the fact that the electrochemical potential re- 
mains constant along the nanotube, i.e. fj,{z) + eip{z) = 0, 
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where (p{z) is the electrostatic potential 

1 rh/2 

^(^z) = -Fz+- dz' p{z')<^{z - z'), (2) 

which is created by external field and by induced charges. 
The kernel, <I>(z — z'), in Eq. ||2Jl takes a simple form in 
the case of isolated NT, lying on substrate with dielectric 
constant, e 



^e* Jo 



da 



+4r2sin^(a/2)] 



1/2 



(3) 




where e* = (e + 1) /2. With the help of Eqs. ^ and 101, 
the condition of constant electrochemical potential can 
now be presented as a closed integral equation for p{z) 



2 „h/2 

+ dz'p{z')<i>{z-z'), (4) 

J-h/2 



where we assumed z to be positive. Eq. should be 
complemented by the obvious condition that p(z) is odd. 

ASYMPTOTIC SOLUTION OF EQ. (H) 

In order to make use of the small parameter r/h, we 
rewrite the integral on the rhs of Eq. Q as follows 

dz'p{z')^z-z')=^ / dz'^^^^K{z,z'), (5) 

-h/2 Jo OZ' 

where the function K{z, z') is defined as 



rh/2 

K{z,z') = / dzi 



$(z-zi) -$(z + zi) . (6) 



Our most important observation is that, in the limit 
h ^ r, the function K(z, z') can be replaced by 
(2e/£*) ln(/i/4r)9(z - z'), where 6(x) is the step- 
function. Possibility of such a replacement is illustrated 
in Fig. 1. A simple from of K{z, z') allows for a drastic 
simplification of Eq. which transforms from inter- 
gral equation to a simple algebraic (quadratic) equation. 
In particular case. Eg = (metallic NT), we obtain the 
following result for the induced charge distribution 



£*gFz 



l + 2g\n{h/r)' 



(7) 



where we have introduced a dimensionless interaction pa- 
rameter g = 2Ne'^/ne*hvo. The above result for p{z) has 
a logarithmic accuracy, in the sense, that numerical fac- 
tor in the argument of logarithm is not specified. In par- 
ticular, the height of the step-function in Fig. 1 contains 
1/4 under the logarithm. Another contribution to the 
argument of logarithm comes from the smearing of the 
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FIG. 1: Plot of the function K{z,z') [scaled with the factor 
2e ln(/i/4r)/e*] for h/r — 10'^ and different values of z: z = 
O.lh (dotted line), z — 0.2h (solid line), z = 0.3/i (dashed 
line), z = OAh (dash-dotted line). 



step- function in Fig. 1. This smearing gives rise to the 
relative correction ~2zd\n p{z)/dz to In^h/Ar); in partic- 
ular case of Eq. |7J) this correction is equal to —2. Overall, 
the condition of applicability of Eq. is ln(/i/r) ^ 1, 
which is met in most of the realistic situations. It follows 
from Eq. {Tj) that the polarizability, X: of the NT, defined 
as P{F) = dzzp{z) ~ has the form 



X 



£*gK^ 



12 l + 2gln(/i/r-) 



(8) 



It is clear from Eq. (jSj) that the product gh\{h/r) w 
(1.74iV/£*) ln(/i/r) is a ciuantitative measure of the 
"metanicity"of the NT. In the limit of a "long" NT, 
when the product g\n{h/r) is large, we have x = 
£*fe^/24 ln(fe/r), which coincides with the textbook ex- 



InOi/ 





pression |1YI | for polarizability of a perfectly conducting 
ellipsoid with axes r and h ^ r. In this limit, with ex- 
ternal field parallel to the NT axis, the resulting field 
at the NT surface is normal to this surface. The oppo- 
site limit, 2gh\{h/r) <C 1, of a "short" nanotube, when 
external field is altered weakly by the induced charges, 
cannot be achieved even for high dielectric constant of 
the substrate, e.g., e* « 6 for Si. 

Consider now a semiconducting (or strained metallic 
[ill) NT with finite Eg. It is seen from Eq. Q that 
charge separation occurs only when the external field is 
strong enough, namely, F > Fth = Eg /eh. It also follows 
from Eq. that, as F increases, electrons and holes 
emerge at the NT ends, while the strip |z| < Eg/2eF in 
the center of NT remains neutral. The behavior of p{z) 
outside the strip is given by 



e^Eg 
2e 



g^ \n\h/r)^ 



■\i-i{h/r) 




(9) 
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and is illustrated in Fig. 2. From the edge of the neutral 
strip to the "bulk" of NT the slope of p{z) decreases by 
a factor 2gln{h/r)/[2gln{h/r) + 1]. Using Eq. ©, one 
can calculate the induced dipole moment, P{F), of semi- 
conducting NT. Obviously, for F ^ Fth it is the same 
XF as for metallic NT. In the vicinity of the threshold, 
{F — Fth) ^ Fth, the induced charge density not only 
occupies small region near the tips, but is also small in 
magnitude. Therefore, P{F) is quadratic in {F — Fth) 
near threshold, namely 

, , e*h^Fth f F V 



of our basic equation Quantum correction to zt 

comes from the penetration of electronic wave functions 
into the classically forbidden region inside the neutral 
strip. Using Zener's formula for tunneling exponent in 
the Dirac spectrum J^l, the relative quantum correction 
can be presented in the form, Sz'^ / zt ~ (Ib /gzt)"^^^, where 
Ib = /Eg£* has the meaning of the exciton Bohr ra- 
dius. Thus, the condition of smallness of the quantum 
correction to zt coincides with the condition that F is 
weaker than the ionization threshold for an exciton. 



F = Q; CHARGED METALLIC NT 
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FIG. 2: Charge density distribution induced by external field, 
F, along the metallic (dashed line) and semiconducting (solid 
line) NT. 



Experimental situations in which electrostatics of a 
charged NT is important are listed in Ref. Clas- 
sical electrostatic analysis of a charged NT was carried 
out in this paper only for a short NT (with h/r w 4). 
Below we find the distribution of charge analytically in 
the limit r. For a charged metallic NT the condition 
of a constant electrochemical potential takes the form 



M = 



2eN 



h/2 



h/2 



dz'p{z')^{z - z'). (12) 



Now p{z) is even. To employ the above ansatz, we take 
derivative from the both sides of Eq. (|12|) and perform 
integration by parts in the rhs. This yields 



SMALL Eg; FINE STRUCTURE OF THE 
NEUTRAL STRIP 

The boundaries of the neutral strip, zt = zkEg/2eF, 
were found from Eq. (0J within the "local" approxima- 
tion. There are two sources of corrections to this re- 
sult: classical and quantum. Classical correction origi- 
nates from the fact that for small Eg the neutral strip is 
surrounded by long charged regions with opposite signs 
of charge. Setting z = zt in Eq. Q), and substituting 
into the rhs the zero-order result © for p{z), we obtain 
the following modified equation for zt 



Zt 



2eF 



2g\n{h/zt) 
1 + 2g\ii{h/r 



■Zt- 



(11) 



The second term in the rhs of Eq. (|ll|l reflects the fact 
that potentials, created by the left and right charged 
neighbors of the neutral strip, do not compensate each 
other completely. As follows from Eq. Hll|l . the rela- 
tive classical correction to the position of the boundary, 
Zt, is bzfjzt w 2 In (/i/zt)/ In (/i/?-). This correction is 
small provided that eFr ^ Eg. The latter condition 
also insures that the underlying energy spectrum of the 
NT is not affected by F ^ as was assumed in derivation 



p(hl2) [$_ - = 



eR{z) 



gs 



h/2 

dz'R{z')^{z-z'), (13) 

-h/2 



where R{z) = dp{z)/dz is the odd function of z, and 
a short-hand notation $±(z) = <I>(-| ± z) is introduced. 
Now Eq. ((T^ has the form similar to Eq. Q . It should be 
complemented by the condition 2 J^^'^ dzp{z) = Q, where 
Q is the total charge on the NT. Using the similarity 
between Eqs. (|13|l and 10}, we readily obtain 



P{z) = ^ 



1 + .g ln(/i/r^ 



■In 



/i2 - 4z2 



(14) 



Hence, p{z) is essentially constant along the NT 0, and 
raises sharply only near the tips. This behavior compares 
favorably with numerical results |l5l | . Logarithmic diver- 
gence in Eq. (|14|l is terminated in the vicinity of the tip 



h/2 



- -y < 



r, so that the net growth of p{z) is given by 



p{h/2)/pi0) = [1 + 2gHh/r)]/[l + gHh/r)]. 

In optical experiments on separated NTs the tubes usu- 
ally have wiggly shapes. Then their response to the exter- 
nal field can be quite peculiar, with numerous alternating 
positively and negatively charged regions separated by 
neutral strips (see Fig. 3a) . An insight into electrostatics 
of a wiggly NT can be drawn from a model example that 
allows for exact solution. 
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BENT NT 

We consider a NT in the form of a semicircle of a ra- 
dius, R. In two limiting cases, the electric field is either 
pointed along the diameter connecting the NT tips (par- 
allel geometry) or perpendicular to this diameter (per- 
pendicular geometry, see Fig. 3). Both geometries are 




FIG. 3: Schematic illustration of charge separation in semi- 
conducting NT in external field: (a) NT of a wiggly shape; 
(b) and (c) NT of a semicircle shape with F parallel and per- 
pendicular to the diameter, respectively. 

described by Eq. Q written in polar coordinates 



eFt{e)R -C^ sign(p)^i?2/4 + [ep{e)/ge*Y 
eR r-/^ 



+ - 



J-Tv/2 



$ i?sin- 



(15) 



where Ft{9) is the tangent component of F. The con- 
stant, C, which determines the electrochemical potential 
of NT, must be found from the condition of the net NT 
neutrality, dOpiO) = 0. We emphasize that in both the 
parallel and perpendicular geometries only the tangent 
component of electric field is responsible for the charge 
separation. The normal component field has no effect 
on the charge separation if the condition Eg 3> eFr is 
satisfied. The latter condition means that field-induced 
mixing of transverse subbands is negligible. 

(i) In parallel geometry, we have Ft = i^sin6'. Thus, 
p{9) is odd, so that C = 0. Same ansatz as for a straight 
NT, yields the following solution of Eq. ((T^ for Eg ~ 
p{d) ^ e*gFRsme/[l + 2.g ln(i?/r)]. It is also easy to 
see from Eq. (|15|) that, for finite Eg, the neutral strip 
occupies the segment 16*1 < arcsin (2eFi?/£'g) of the NT 
near its top (Fig. 3b). 

(ii) In perpendicular geometry. Ft (9) = Fcosd, is the 
even function of 9. As a result, p{9) differs qualitatively 
from the case of a straight NT. The form of p{9) for this 
geometry can be found from the same ansatz Eq. H13|) 
that was used for the straight charged NT. In particular, 
for £'g = we obtain 



p{9) = e*FR 



g (cos 61 - 2/7r) 
l + 2.gln (R/r)' 



(16) 



The two points at which p{9) changes the sign are there- 
fore located at 6* = ±6*0 = ±50.5°. For a finite Eg two 



neutral strips are formed around 9 = ±9q. Their bound- 
aries, ±0+ and ±0_, (see Fig. 3c) arc determined by the 
conditions 



2C + E, 



a 



2eFR 



cos 9^ 



2C-E, 



a 



2eFR 



(17) 



When the gap is small. Eg -C eFR, the centers of the 
strips arc still located near ±9o, while the strip width, 
69 ^ {9_ - 9 + ), can be found directly from Eq. H17|l . 
namely 69 = Eg/eFRsin9o « l.iEg/eFR. 

For small ratio Eg/2eFR it is also easy to trace 
the crossover between the parallel and perpendicular 
geometries as the field is rotated. For rotation angle, 
/?, and for Eg = 0, a, straightforward generalization 
of Eq. (|16|l yields the following angular dependence 
of the charge density p{9) ex [cos(/3 — 9) — ^ cos /?] . 
From this dependence we conclude that, as f3 decreases 
from 90° (parallel geometry), the narrow neutral strip 
at the top of semicircle moves to the left. At critical 
Pc = arctan(2/7r) w 32.5°, when the neutral strip is 
located around 9 w 25°, the second neutral strip emerges 
at the right end of the NT. As (3 decreases further, 
both neutral strips move to the left and assume their 
"perpendicular" positions 9 = ±50.5°. 



EXPERIMENTAL IMPLICATIONS 

In Ref. I23 the potential distribution along the NT 
was measured using the atomic force microscope. Ex- 
perimentally measured profile of the voltage drop in the 
structures with a small contact resistance is consistent 
with existence of a neutral stri p n ear the NT center. 
Note, that the substrate in Ref. |20| was thick: D = 200 
or D = 500nm, but still thinner than the NT length, 
h — 1200nm. Taking into account the presence of a gate 
at distance D 5 , h from the NT amounts to replacement 
m m m m of \n{h/r) by \n{D/r) in all the above 
formulas. 

Our results have direct relevance to the electro-optics 
of NTs. Measurements of electroabsorption in single- 
walled carbon NTs were recently reported in Ref. [la- It 
might seem that with photon energy ~ leV much bigger 
than Eg the large-scale nanotube geometry is not impor- 
tant. This, however, is not the case. The reason is that 
the dipole moment of the many-body optical transition 
[il,!!^ is directed along the tube 0113. Even if exter- 
nal field is parallel to the tube axis, the resulting field is 
almost perpendicular to the nanotube surface, and thus 
its effect on the optical transitions is suppressed. This 
strong suppression by a factor [1 + 2ghi(h/r)] must be 
taken into account when the oscillator strength is ex- 
tracted from electroabsorption On the other hand, 
within the neutral strip, the acting field is equal to the ap- 
plied field. However, the relative contribution, Eg/eFh, 
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of neutral strips to the absorption is small. 
As a final remark, note that dimensionless parameter, g, 
which governs the screening properties of the NT, has a 
transparent meaning. For 2D electron gas with a den- 
sity of states, i^, the linear screening length is equal to 
I = e* I2-Kve^ . If the gas is wrapped into a cyhnder of a 
radius, r, then the degree of penetration of external field 
inside the cylinder is determined by the ratio rjl. Up 
to a numerical coefficient, this ratio is nothing but the 
parameter g. 
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